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Abstract. For the 3-component dispersionless Boussinesq-type system, we construct 
two compatible nontrivial finite deformations for the Lie algebra structure in the sym- 
metry algebra. 



Introduction. In this short note we construct a two-parametric family of nontrivial 
finite deformations for the Lie bracket in the algebra of symmetries for the 3-component 
dispersionless Boussinesq system of hydrodynamic type [U [2] 

^ ■ S = {u t = ww x + v x , v t = -uw x - 3u x w, w t = u x | . (1) 

First, we establish a nontrivial property of a previously known [2] self-adjoint Noether 
operator Aq\ cosym£ — > sym£: its image is closed w.r.t. the commutation. Hence this 
operator, and the bi-Hamiltonian pair (Ai,A 2 ) for (CD), see [lj, transfer the standard 
bracket [ , ] in sym S to the Lie algebra structures on their domain. We prove that the 
three new brackets are compatible. 

The Noether operator A is invertible on an open dense subset of S. This yields 
^ • two recursion operators Ri = Aj o Aq . sym£ — > sym£ . The images of Ri are again 

closed w.r.t. the commutation, and this property is retained by their arbitrary linear 
! combinations. Using the 'chain rule' formula ( TTT1) for the bi-differential brackets on 

domains of the operators A4 and Ri, we calculate the second Lie algebra structures 
[, } Rt on sym£. 

All notions and constructions are standard [31 Hj. We stress that the concept of line- 
ar compatible differential operators with involutive images, which we develop here, can 
be applied to the study of other integrable systems with or without dispersion (e.g., 
see 0). 

Remark 1. Let M be a smooth finite-dimensional orientable real manifold. The con- 
struction of trivial infinitesimal deformations 

[x,y]x ■= ^| A=0 e " AN K N (x),e AN (y)] = [Nx,y] + [x,Ny] - N([x,y}) 

of the standard Lie algebra structure [ , ] on the tangent bundle TM 3 x, y is well- 
developed in the literature [6]: if the Nijenhuis torsion \Nx, Ny] — N([x,?/]n) for an 
endomorphism N: T(TM) — > T{TM) vanishes, then the Lie brackets [, obtained 
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by iterations N fc of the Nijenhuis recursion N are compatible (their linear combinations 
are Lie algebra structures as well). Suppose further that M is equipped with a Poisson 
bi- vector V G T(/\ 2 (TM)). If the Nijenhuis and Poisson structures (N, V) satisfy 
two compatibility conditions [B], then they generate infinite hierarchies of compatible 
Poisson structures N fc o V, k > 0. The concept of Poisson-Nijenhuis structures admits 
a straightforward generalization [7] for the infinite jet bundles over smooth manifolds 
and for infinite-dimensional integrable systems of PDE. 

On the contrary, in this paper we construct two nontrivial finite deformations [, 
of the standard Lie bracket [ , ] on the symmetry algebra sym£ for ([T]). We shall use two 
local recursion operators Ri, i = 1,2, whose images are closed w.r.t. the commutation. 
Thence we obtain the new bracket through 

[Ri^Px, Riip 2 ] = Ri([fi, ^ 2 ]r z ) for any (fi,(f 2 e sym£. (2) 
The construction constits of two steps. 

Involutive distributions and linear compatibility. First, let a linear operator □ 
in total derivatives be either a recursion sym £ — > sym £ for an evolutionary system £, 
or a Noether operator cosym £ — > sym £ whose arguments are the variational covectorgj 
for £ = {F = 0} ('cosymmetries' ip = V*(l) such that divt] = (1, V(F)) = (ip, F) for a 
conserved current r\ = pdx+- ■ ■ and the adjoint V* of the arising operator V, see pill]). 
For example, all Hamiltonian operators for £ are Noether. 

Suppose further that the image of □ is closed w.r.t. the commutation in sym£: 
[imD,imD] C imQ The Lie algebra structure [, ]|. . is transferred by □ onto the 
quotient Q = domD/kerD: 

[□(00, □(0")] =n([0',0 /, ] n ), ^/en. 

By the Leibnitz rule, two pairs of summands appear in the bracket of the evolutionary 
vector fields fields £.n(<j>>) and £□(</,»): 

[□(0'), □(0")] = □(£□(*) - £□(*») + (£n(^)( D )(^) - £n(0")( n )(0O)- 
In the first summand we have used the permutability of evolutionary derivations and 
operators in total derivatives. The second summand hits the image of □ by construction. 
Therefore the bracket [0', 0"]n equals 

[(f)', 0"]n = £□(<//) (0") - £□(<*.») (00 + {{0'> 0"}}n- ( 3 ) 
It contains the two standard summands and the skew-symmetric bilinear bracket {{ , }}□. 

Example 1. A self-adjoint (hence non-Hamiltonian) zero-order Noether operator Aq 
for ([[]) was found in [2j: 

ww x + v x —3u x w — uw x 

Ao = | — 3u x w — uw x —3w 2 w x — Av x w — uu x v x ] . (4) 



We label the evolution equations upon u 1 , . . . , u m in the system E with the same variables u % that 
occur in the left-hand sides. With such convention, the cosymmetries "0 equal the variational derivatives 
5p/Su of the conserved densities p. This yields the transformation law for Noether operators under 
reparametrizations of u l that preserve the evolutionary form of £ . 
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We discover that its image in sym£ is involutive; the components of the arising bracket 
{{p, q}}Ao with p , q e cosym£ are 

{{P, q}} U A = V W X Q U - P U Q W X + MP U Q V X ~ f x q u ) + Mp v ti ~ vl<?) 

+ P u x q w ~ P w q u x + 2w x (p u q v - p v q u ) + u(p v q v x - plq v ), 

m q}} v Ao = p u x q u - p u q u x + Mp v q v x - pW) + p v x q w - p w ql + p w x q v - p v q w x , 

{{P, q}}A = u{p u q v x - p v x q u ) + Zw 2 (fql - p v x q v ) + 2u x (fq u - p u q*) 
+ w{p u x q u - p u q u x ) + u{p v q u x - p u x q v ) + p w x q w - p w q™. 

Lemma 1 ([HI H]). The image of any Hamiltonian operator A = ||^ r ^4^ ■ -D r || is 
closed w.r.t. the commutation; the bracket {{ , }}^ on its domain equals 

{{«W}i = £(-i)'(^° [E^W-^])(«- (5) 

Example 2. The pair [A\,A-i\ of compatible Hamiltonian operators for (0Q) was ob- 
tained in pQ: 

/ D x \ 

A x = -4wD x - 2w x D x , (6) 
V D x / 

'{2w 2 + Av) D x + 2{ww x + v x ) —Wuw D x — (5u x w + 9uw x ) 3u D x + u x 
A 2 = | —lluwD x — Qu x w — 2uw x 2hD x + h x Av D x + v x 

3u D x + 2u x Av D x + 3v x 2w D x + w x 

' (7) 

where we put h = — (|u 2 + Svw + 3w 3 ). The components of the brackets {{ , }}^. are 
given by ([Sj): for any p , q e cosym£ they equal, respectively, 

{{ , }}\ = {{ , }}\ = 0, {{P, q}}\ = 2 {p v q v x - PW) (8) 

and 

m q}}% = 2 (K<f ~ p u £) + MpX - pIi u ) + ^(p v q: - ptq v ) 

+ p u x q w - p w q u x + Aw x (p u q v - p v q u ) + 3u{p v q v x - p v x q v ), (9a) 

{{p, q}}\ = mq u - p u q u x ) + M« - pW) + f x q w - P w ql + Kp w x q v - p v q w x ), (9b) 

{{P, q}}'l 2 = Mp U ql - P V S U ) + (8w + 9w 2 ) ■ {fql - p v x q v ) + Au x {fq u - p u q") 

+ 2w(p u x q u - p u ql) + 9u(p v q u x - p u x q v ) + p w x q w - p w q w x . (9c) 

We claim that the Noether operator Aq, whose image in sym£ is closed w.r.t. the 
commutation, is compatible with A\ and A 2 in this sense. 

Definition. We say that N > 2 operators A^: (co)sym£ — > sym£ on £ with a com- 
mon domain and involutive images, [im Aj, im Aj] C imAj for 1 < i < N, are linear 
compatible if their linear combinations A^ = Yli=x retain the same property of 
involutivity for any A. 
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Example 3. It can easily be checked that three Noether operators (j3J), ©, and ([7j) for 
system ([T]) are linear compatible. 

We note that linear compatible Hamiltonian operators are Poisson compatible, and 
vice versa, because formula (jSJ) is linear in coefficients of A. 

Theorem 2. The bracket {{ , }}a^ on the domain of the combination of linear 
compatible operators A{ is 

N 

i=l 1=1 

The pairwise linear compatibility implies the collective linear compatibility of Ax, . . . , A^. 

Proof. This is readily seen by inspecting the coefficients of A? in the quadratic polyno- 
mials in \i that appear in both sides of the equality [A^(p), ^x(^)] = ^a([P' Q] A x)' ^ ere 
p,q G f2 = dom A*/ fljLi ker Aj for any i. 

Consider the commutator AiA(p), ^ - XjAj(q)j . On one hand, it is equal to 

= X>Ai[A(p), MQ)] +J2 X2 M^A i ( P )(q) ~ Za M (p) + {{P,Q}}A). (10) 

ijtj i 

On the other hand, the linear compatibility of Ai implies 

= A X (E Ax(p) (q)) - A x (8. Axiq) (p)) +A x ({{p,q}} Ax ). 

The entire commutator is quadratic homogeneous in A, whence the bracket {{ , }}a^ is 

linear in A. From (flOl we see that the individual brackets {{ , j}^ are contained in it. 
Therefore, 

{{p, q}}a x = ^2^e- {{P, q}}a ( + E Xe ' 7£ ( p ' q ^ 
i t 

where 7^ : Q x Q — > fi. We claim that all summands 7^(-, •), which do not depend on 
A at all, vanish. Indeed, assume the converse. Let there be £ G [1, . . . , N] such that 
liiPiQ) 7^ 0; without loss of generality, suppose I = 1. Then set A = (1,0,..., 0), 
whence 

[E A ^(p)>Z)'W?) 

+ (A1A1) ^£(AiAi)(p)(?) - £(AiAi)(ff)(p) + M{P:?}UiJ- 

Consequently, 7^(0, g) G ker for all p and g. Now we are forced to use the nondegen- 
eracy assumptions P| T kerA£ = {0} for the operators A e = || ^ r A\ ■ D T \\. Finally, we 
have that 7^ = for all £, which concludes the proof. □ 



^We assume that there is no functional freedom in the kernels of Ag. This restriction and its 
consequences will be discussed in full detail in a separate paper. 



(X 1 A 1 )(p),(X 1 A 1 )( q ) = (A 1 A 1 )(A l7 i(p,g)) 
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Corollary 3. Two such operators A and B with involutive images are linear compatible 
iff for any p, q G f2 one has 

[B(p),A{q)] + [A{q),B(p)] = A(\p, q] B ) + B(\p, q) A ), 

which is equivalent to the relation 

E A{p )(B)(q) + 8, B(p) (A)(q) - E Aiq) (B)(p) ~ ZB {q) {A){p) = A({{p,q}} B ) + B ({{p, q}} A ) . 

Recursion operators and new Lie brackets. The second step in the construction 
of the new Lie brackets on sym£ is as follows. We note that the zero-order opera- 
tor (jlj) has the inverse lu = Aq 1 on an open dense subset of £. This yields the local 
recursion operators Ri := Aj o A : sym£ — > sym£. By construction, their images are 
closed w.r.t. the commutation [, ]. (We note also that the images of the Hamiltonian 
operators A^ for S are generally not closed w.r.t. [, here 1 < i, j < 2.) Therefore 
we determine the new Lie algebra structures [, ]n i on sym£ by (j2J). Let us summarize 
the result. 

Proposition 4. The recursion operators Ri = A\ o A^ 1 , i = 1,2, for system (PQ) 
are linear compatible. The new Lie algebra structures [, } Bi on sym£ span the two- 
dimensional space of compatible nontrivial finite deformations of the standard bracket [ , ] 

The transformation rules for Ri and hence for [, under any reparametrizations of 
the variables u,v,w, are obvious. The arising bi-differential brackets {{ , }}^, see (EJ), 
are obtained from {{ , }} A , using the following 'chain rule' (c.f. [8]). 

Theorem 5. Suppose that the images of linear differential operators A and R = A o u 
are closed w.r.t. the commutation of evolutionary vector fields. Then the brackets {{ , }} A 
and {{ , }} R are related by the formula 

= £h(6)H(6) - £jk6)(w)(£i) + (11) 

for any sections £1,^2 that belong to the domain of u. 

Proof. Denote ipi = k>(&) and ipi = A(ipi) for i = 1,2. We have 

[ip ll<P2 ] = (Ao W )(6^fe) - 8^(6) + (12a) 
On the other hand, we recall that ^ = k>(£i) and deduce 

= a(£^(^ 2 ) - e V2 (Vi) + M^2}U) (12b) 

= (io W )(£^fe) - £ V2 (6)) + A(£ V1 H(6) - £ w (o;)(ei) + {{ikMU). 
Now subtract (11 2 aft from (112bl) . Omitting the operator A, we obtain the assertion. □ 

The bi-differential brackets {{ , }}r { for the recursions Ri constructed above are com- 
pletely determined by the chain rule (fTTj) with uj = Aq 1 : 

{{fi,?*}}* = Ao(£fli( w )(w)(^ 2 ) - £^( V2) (o;)(^i) + {{wW.wW}^), (13) 

here (pi, (p% G sym£ are any symmetries of ([T]). The three components of each bracket 
{{ , }}i^ can be calculated explicitly, e.g., using the environment [9], see Appendix lAl 
The coefficients of the skew-symmetric couplings D%(ipi) ■ D^(<f\) — D%((p\) ■ D^ip^), 
< a + (3 < 1, 1 < a, b < 3, are relatively big due to the presence of the powers 
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(detA ) a , 1 < a < 3 in the denominators. Be that as it may, the two local recursion 
operators Ri generate the first known examples of compatible well-defined new Lie 
brackets on the symmetry algebra of system ([T]) via the Yang-Baxter equation 

Remark 2. The application of the classical r-matrix formalism pU] for a given Lie 
algebra g generates Liouville integrable systems using the second Lie algebra structure 
[, ] r that solves the Yang-Baxter equation YB(a). Since equation ([2]) corresponds to 
the degenerate case a = 0, we pose the problem of finding recursion operators for 
= sym 8 that will lead to relevant factorizations and produce new integrable systems. 

Remark 3. Particular examples of non-Hamiltonian linear operators with involutive im- 
ages, and Lie brackets on their domains, are scattered in the literature (e.g., see [TT] 
for a dispersionless set-up and a bracket of 1-forms). However, let us remember that, 
first, the differential order of such operators can be sufficiently high for systems with 
dispersion. Second, the domains and images of such operators can be composed by 
(co) symmetries of two different equations. For example, a class of higher-order oper- 
ators with involutive images is known for the open 2D Toda chains and the related 
KdV-type systems, see [5]. Involutive distributions of operator- valued evolutionary 
vector fields will be the object of a subsequent publication. 
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Appendix A. Calculation of {{ , }} Ri for the two recursions Ri 

The following program for the Jets environment [H] under Maple calculates the 
bracket {{ , }} Rl on the domain of the recursion operator R 1 — A\ o Aq 1 , see @ and (jSJ), 
for the dispersionless 3-component Boussinesq-type system ([T]). The bracket {{ , }}^ 2 
induced by R2 is obtained using a slight modification of this program. The modification 
amounts to a substitution of A 2 and {{ , }}_| , see (jHJ), for A\ and {{ , }}_| , respectively. 

> read ' Jets . s 1 ; 

> coordinates ( [x] , [u, v, w,pl ,p2,p3,ql , q2 ,q3] ,5) : 

The dependent coordinates p l , q 7 denote the components of the sections <fi,<f2 in the 
domain sym£ of R\. The notation is in parallel with formula ( TT31) . 

> A : =Matrix (3,3): ia : =Matr ix (3,3): M : =Matrix (3 , 3) : N : =Matrix (3,3): 

> A : =<<w*w_x+v_x , -3*w*u_x-u*w_x , u_x> I 

<-3*w*u_x-u*w_x , -3*w"2*w_x-4*w*v_x-u*u_x , v_x> I <u_x , v_x , w_x>> ; 

We have assigned A = A , see (@|. 

> with(LinearAlgebra) : 

> ia:=MatrixInverse(A) : 

This yields ia = Aq 1 : sym£ — > cosym^, which we also denote by u>. Now put m = 
Aq 1 ^) and n = Aq 1 ^). 

> ml := simplify(ia[l,l]*pl+ia[l,2]*p2+ia[l,3]*p3) : 

> m2 := simplify(ia[2,l]*pl+ia[2,2]*p2+ia[2,3]*p3) : 

> m3 := simplify(ia[3,l]*pl+ia[3,2]*p2+ia[3,3]*p3) : 
> 

> nl := simplify(ia[l,l]*ql+ia[l,2]*q2+ia[l,3]*q3) : 

> n2 := simplify(ia[2,l]*ql+ia[2,2]*q2+ia[2,3]*q3) : 

> n3 := simplify(ia[3,l]*ql+ia[3,2]*q2+ia[3,3]*q3) : 

We have computed ipi — ^(fi)', we recall that i[) G cosym£. 

Next, we apply the Hamiltonian operator A\ and obtain symmetries of (PQ), k = 
ii(m) = A 1 (Ao l M) an d 1 = ii(n) = A 1 (A^ 1 (<p 2 )) . 

> kl := simplify(evalTD(TD(ml,x))) : 

> k2 := simplify(evalTD(-4*w*TD(m2,x)-2*w_x*m2+TD(m3,x))) : 

> k3 := simplify (evalTD(TD(m2,x))) : 
> 

> ll:=simplify(evalTD(TD(nl,x))) : 

> 12:=simplify(evalTD(-4*w*TD(n2,x)-2*w_x*n2+TD(n3,x))) : 

> 13:=simplify(evalTD(TD(n2,x))) : 
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Now we act by evolutionary derivations on the coefficients of the operator u = A Q , 
see (ITT]) and ( fl3l) . We set M = ^r x {u, x ) (ijS)\ note that the matrix M is symmetric. 

> M[l,l] :=simplify(evalTD(kl*pd(ia[l,l] ,u)+TD(kl,x)*pd(ia[l,l] ,u_x)+ 
k2*pd(ia[l,l] , v)+TD(k2,x) *pd(ia[l , 1] ,v_x)+k3*pd(ia[l , 1] ,w)+ 
TD(k3,x)*pd(ia[l,l] ,w_x))) : 

> M[l,2] :=simplify(evalTD(kl*pd(ia[l,2] ,u)+TD(kl ,x)*pd(ia[l ,2] ,u_x)+ 
k2*pd(ia[l,2] ,v)+TD(k2,x)*pd(ia[l ,2] ,v_x)+k3*pd(ia[l ,2] ,w)+ 
TD(k3,x)*pd(ia[l,2] ,w_x))) : 

> M[l,3] :=simplify(evalTD(kl*pd(ia[l,3] ,u)+TD(kl ,x)*pd(ia[l ,3] ,u_x)+ 
k2*pd(ia[l,3] ,v)+TD(k2,x)*pd(ia[l,3] , v_x)+k3*pd(ia[l ,3] ,w)+ 
TD(k3,x)*pd(ia[l,3] ,w_x))) : 

> M[2,2] :=simplify(evalTD(kl*pd(ia[2,2] ,u)+TD(kl,x)*pd(ia[2,2] ,u_x)+ 
k2*pd(ia[2,2] ,v)+TD(k2,x)*pd(ia[2,2] ,v_x)+k3*pd(ia[2 ,2] ,w)+ 
TD(k3,x)*pd(ia[2,2] ,w_x))) : 

> M[2,3] :=simplify(evalTD(kl*pd(ia[2,3] ,u)+TD(kl ,x)*pd(ia[2,3] ,u_x)+ 
k2*pd(ia[2,3] , v)+TD(k2,x) *pd(ia[2,3] ,v_x)+k3*pd(ia[2,3] ,w)+ 
TD(k3,x)*pd(ia[2,3] ,w_x))) : 

> M[3,3] :=simplify(evalTD(kl*pd(ia[3,3] ,u)+TD(kl ,x)*pd(ia[3,3] ,u_x)+ 
k2*pd(ia[3,3] , v)+TD(k2,x) *pd(ia[3,3] ,v_x)+k3*pd(ia[3,3] ,w)+ 
TD(k3,x)*pd(ia[3,3] ,w_x))) : 

In the same way, we define the symmetric matrix N = E r^^uj) . 

> N[l,l] :=simplify(evalTD(ll*pd(ia[l,l] ,u)+TD(ll,x)*pd(ia[l,l] ,u_x)+ 
12*pd(ia[l,l] ,v)+TD(12,x)*pd(ia[l,l] ,v_x)+13*pd(ia[l,l] ,w)+ 
TD(13,x)*pd(ia[l,l] ,w_x))) : 

> N[l,2] :=simplify(evalTD(ll*pd(ia[l,2] ,u)+TD(ll ,x)*pd(ia[l ,2] ,u_x)+ 
12*pd(ia[l,2] ,v)+TD(12,x)*pd(ia[l,2] ,v_x)+13*pd(ia[l ,2] ,w)+ 
TD(13,x)*pd(ia[l,2] ,w_x))) : 

> N[l,3] :=simplify(evalTD(ll*pd(ia[l,3] ,u)+TD(ll ,x)*pd(ia[l ,3] ,u_x)+ 
12*pd(ia[l,3] ,v)+TD(12,x)*pd(ia[l,3] ,v_x)+13*pd(ia[l ,3] ,w)+ 
TD(13,x)*pd(ia[l,3] ,w_x))) : 

> N[2,2] :=simplify(evalTD(ll*pd(ia[2,2] ,u)+TD(ll,x)*pd(ia[2,2] ,u_x)+ 
12*pd(ia[2,2] ,v)+TD(12,x)*pd(ia[2,2] ,v_x)+13*pd(ia[2 ,2] ,w)+ 
TD(13,x)*pd(ia[2,2] ,w_x))) : 

> N[2,3] :=simplify(evalTD(ll*pd(ia[2,3] ,u)+TD(ll,x)*pd(ia[2,3] ,u_x)+ 
12*pd(ia[2,3] , v)+TD(12,x) *pd(ia[2,3] ,v_x)+13*pd(ia[2,3] ,w)+ 
TD(13,x)*pd(ia[2,3] ,w_x))) : 

> N[3,3] :=simplify(evalTD(ll*pd(ia[3,3] > u)+TD(ll,x)*pd(ia[3,3] ,u_x)+ 
12*pd(ia[3,3] , v)+TD(12,x) *pd(ia[3,3] ,v_x)+13*pd(ia[3 ,3] ,w)+ 
TD(13,x)*pd(ia[3,3] ,w_x))) : 

We act by the operators M, N on (f2 and ipi, respectively, and calculate the difference 
e = filial) MM - £ri(<p 2 )(u)M- 

> el : =simplif y (M [1 , 1] *ql+M [1 , 2] *q2+M [1 , 3] *q3-N [1,1] *pl-N [1 , 2] *p2-N [1,3] *p3) 

> e2 : =simplif y (M [1 , 2] *ql+M [2 , 2] *q2+M [2 , 3] *q3-N [1,2] *pl-N [2 , 2] *p2-N [2 , 3] *p3) 

> e3 : =simplif y (M [1 , 3] *ql+M [2 , 3] *q2+M [3 , 3] *q3-N [1 , 3] *pl-N [2 , 3] *p2-N [3 , 3] *p3) 
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Next, we recall that the bracket {{ , }}^ for A\ is equal to (jHJ), and substitute = 
co((pi) in it. Thus we put s = e + {{m, n}}* . 

> si :=simplify(el) : 

> s2:=simplify(e2) : 

> s3:=simplify(e3+evalTD(2*(m2*TD(n2,x)-TD(m2,x)*n2))) : 

We can now check that formulas (jSj) is correct. 

> J:=Jacobi([kl,k2,k3,0,0,0,0,0,0] , [11,12,13,0,0,0,0,0,0]) : 

> Dl:=simplify(evalTD(TD(sl,x))) : 

> D2:=simplify(evalTD(-4*w*TD(s2,x)-2*w_x*s2+TD(s3,x))) : 

> D3:=simplify(evalTD(TD(s2,x))) : 

> simplify(evalTD(J[l]-Dl)) ; 



> simplify(evalTD(J[2]-D2)) ; 



> simplify(evalTD(J[3]-D3)) ; 



Finally, we act onto the right-hand side of (jTTI) by the operator uo^ 1 = A Q , and thus 
we obtain the components Z = A (s) of the bracket {{ , j}^. 

> Zl :=simplif y ( (w*w_x+v_x) *sl+(-3*w*u_x-u*w_x) *s2+u_x*s3) ; 

> Z2 : =simplif y ( (-3*w*u_x-u*w_x) *sl+ (-3*w*w*w_x-4*w*v_x-u*u_x) *s2+v_x*s3) ; 

> Z3:=simplify(u_x*sl+v_x*s2+w_x*s3) ; 

The result is somewhat a surprise: the output contains more than 15,000 lines. How- 
ever, decomposition ffl3l) makes the expressions manageable. Obviously, the number of 
essential summands is reduced twice by the skew-symmetry. Next, the summands in the 
three components of {{ , j}^ are collected at the nine bilinear expressions p 1 ■ q^ —p> x ■ q\ 
1 < 2, j < 3 and at the six multiples p l ■ qi —p>-q l with 1 < i ^ j < 3. All the coefficients 
are fractions of differential polynomials in u,v,w with (det^o) 3 in the denominators; 
many of these fractions are reducible. 



